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O . The three neutrino Zee model and its extension includine; three active and one sterile species are 
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Q I studied in the light of new neutrino oscillation data. We obtain analytical relations for the mixing 
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angle in solar oscillations in terms of neutrino mass squared differences. For the four neutrino 



^ ■ case, we obtain the result sin^200 1 — [(Am^^^)^/(4 ArriLSMDAmQ)]^, which can accommodate 

> , 

, both the large and small mixing scenarios. We show that within this framework, while both 

CN ] the SMA-MSW and the LMA-MSW solutions can easily be accommodated, it would be difficult 

to reconcile the LOW-QVO solutions. We also comment on the active-sterile admixture within 
phenomenologically viable textures. 
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I. INTRODUCTION 



Recent results from SNO |l| and super-Kamiokande experiments have confirmed the 
presence of a non-electron flavor in the measured solar v^, flux on earth implying the ex- 
istence of neutrino oscillations. A detailed analysis of the data favors the Large Mix- 
ing Angle (LMA) solution [0] within the Mikheev-Smirnov-Wolfenstien (MSW) framework, 
though other solutions are still not ruled out. A simplified two-flavor oscillation picture 
prefers a neutrino mass squared difference AitIq ~ 4 x 10~^ eV^ with sin^26'0 ~ 0.66 for 
the mixing angle Qq. On the other hand, results Q from the super-K atmospheric neu- 
trino experiments indicate a non-muonic component in the atmospheric flux. Again, in 
terms of the difference of squared masses between two oscillating neutrinos, these imply 
^■^Atm ~ 3 X 10^^ eV^ with a near maximal mixing sin^2^^Atm ~ 1- These two indepen- 
dent scales and mixing angles can be accommodated within the standard model of three 
neutrinos z/g, z/^, z/^. However, in addition to these, another scale pertaining to neutrinos 
is indicated by ^ z7e (and z/^j ^ flavor oscillations seen by the LSND experiment [^. 
The parameters of a two-flavor oscillation that explain these results are ArriLSMD ~ ^(1 ^V^) 
and sin^26'LSND ~ 0(10^^). A simultaneous explanation of all the three types of oscillations 
is not possible with only three types of neutrinos. 

An additional light neutrino, which is not electroweak active and is hence called the sterile 
neutrino z^s, is often proposed to understand all the above anomalies simultaneously. Within 
this four neutrino framework, only two distinct mass patterns (the 3 -|- 1 and the 2 -|- 2) 
were originally allowed by the data. The former has three nearly degenerate active neutrinos, 
with mass differences equal to the solar and atmospheric oscillation scales, separated as a 
whole in mass from the sterile flavor by the large LSND mass difference. However, it is 
not favored when fltted to the latest LSND results, together with data from reactor 
experiments like BUGEY; therefore, we discard it. In the 2 -|- 2 case, two doublets, each 
nearly degenerate with mass splits 0(yAm|) and m^^^) and corresponding mixing 

angles Qq and ^^Atm respectively, are separated by a large mass difference C(^ Aitilsnd). The 
mixing between the two doublets is controlled by the angle 6'lsnd- 

It is true that, strictly within two flavor oscillations, neither the SNO nor the super-K 
results support the occurrence of a sterile species in either solar or atmospheric neutrino 
oscillation. Nevertheless, this simple picture changes when one considers all four neutrinos 
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FIG. 1: The 2 + 2 neutrino mass pattern. The heavier (lighter) pair is spht by the mass scale of 
the atmospheric (solar) neutrino anomaly while the two pairs are separated by the LSND scale. 

Ue, I'fj., i^Ti T^s- A comprehensive analysis within a four neutrino framework of all (including 
the latest) data on solar, atmospheric and LSND oscillations, also taking into account extant 
reactor and accelerator constraints, has recently been performed by Gonzalez-Garcia et.al 
0. To summarize their results, it is convenient to define the four species neutrino mixing 
matrix as 

^ Uel Ue2 Ue3 Ue4 ^ 
U^l U^2 U^3 U^4 
Url Ut-2 Ut-3 Ut-4 

\ u,i u,2 u,3 u,4 ; 

In eq.(|I|), the subscripts (1,2) represent the lighter pair participating in solar neutrino os- 
cillations, while (3, 4) refer to the heavier pair relevant to atmospheric neutrino oscillations. 
These pairs are separated in mass by the LSND mass scale. We present a schematic dia- 
gram of this scenario in Fig.l. For that situation, the BUGEY experiment provides the 
maximum constraint on the content in the heavier (3, 4) pair namely 
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In addition, the CCFR |^ and the CDHSW [IU| experiments constrain the content in the 
lighter (1,2) pair by 
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|U^ir + |U^2n 0.2, (3) 

again for an LSND scale mass separation between the pairs. 

The rest of the constraints on U arrive from solar and atmospheric neutrino oscillation 
data. Eqs. (0) and (|^), show that in the 2 + 2 mass pattern, is mostly confined to the 
(1, 2) pair whereas largely resides in the (3, 4) pair. Three alternative situations are now 
possible |1TI[| : 



(i) Solar neutrino oscillations take place into the purely active neutrino Ur] atmospheric 
neutrinos oscillate into the purely sterile neutrino with maximal mixing. 

(ii) Solar neutrino oscillations take place into the purely sterile neutrino atmospheric 
neutrinos get converted to the purely active neutrino with maximal mixing. 

(iii) Both solar and atmospheric neutrinos oscillate into linear combinations of Us and z/,-; 
the combination pertaining to atmospheric neutrino oscillations is maximally mixed with 

As mentioned earlier, recent data from both the SNO and the super-K ^ experiments 
disfavor two flavor neutrino oscillations into only a sterile species as an explanation for either 
the solar or the atmospheric anomaly. Consequently, configurations (i) and (ii) are severely 
constrained. In contrast, configuration (iii) can still be realized even with new data from 
SNO The favorite 2 + 2 configuration, most favored by the data, is when the linear 
combination of Us and Ur in the state to which z/e from the Sun oscillates, is not maximal 
but in the ratio 1 : 2, i.e, 20 : 80 in the probabilities. We shall define this 'active-sterile 
admixture' as the sterile content in the solar sector: 



A= IU31I2 + IU32I'. (4) 

When the solar mixing angle is within the LMA region, corresponding to the best-fit solution, 
the active-sterile admixture is required to be 0: 

A ^0.18 -0.2. (5) 
This means that the 'atmospheric' pair (3,4) has a dominant sterile content. To be precise. 
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the sterile species participates in about 20% (80%) of the solar (atmospheric) neutrino 
oscillations. It should be emphasized that the SNO results have played a key role in these 
conclusions^ reached by the authors of Ref. |]^. This configuration leads to what has been 
called 10] dose to active solar plus close to sterile atmospheric (CAS + CSA) neutrino 
oscillations. 

In addition to the above 'best-fit' solution, configuration (ii) with the solar neutrino 
oscillations occurring into purely sterile neutrinos is still allowed by the data. The active- 
sterile admixture A in this case is : 

A= |Usi|^ + |Us2|^~0.91-0.97. (6) 

However, here the solar neutrino mixing angle lies in the SMA region. For both of these fits, 
it is seen that 



|U;.i|' + |U^2|'^0. (7) 

In this work, we focus on these mass patterns and study the viability of the 4x4 neutrino 
mass matrix within a radiative Zee-type model |M extended to include the sterile neutrino 



I5|. An extension from the standard 3x3 mass matrix to the 4x4 mass matrix can be 



realized either by a conventional seesaw type of mechanism ||TB[ with heavy right-handed 



states or in a radiative model The three neutrino Zee model has been quite popular 



in analyzing neutrino oscillation data (minus those of the LSND experiment) on account of 
its predictivity. However, this model has run into a serious problem with the data on the 
mixing angle pertinent to solar neutrino oscillations. It has been shown Jl^ that sin'^2^Q 
in this model is forced to be close to unity |18| within (9[(Am|/Am|)^] in disagreement 



with the best-fit value sin^2^^0 ~ 0.66. More specifically, if sin^26'Atm ~ 1 is used as an 



input, the Zee model is found to allow [17, Iql only solutions with bimaximal mixing [19 



This result of the Zee model is a natural consequence of the structure of its mass matrix - 
specifically, its vanishing diagonal elements, rather than the details of the model. It is thus 



^ The active (sterile) content in the solar neutrino sector can be decreased (increased) to more than 
what is implied by eq.(||) if the flux is suitably renormalized. Such a renormalization of up to 30% 
may be allowed by the theoretical uncertainties in the calculation of the said flux. It has been shown 
recently j^, however, that the inclusion of the SK data on the day and night spectral energy distribution 
disfavor such a scenario at 3cr level. 
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essential to understand the consequences of such a mass matrix when extended to the four 
neutrino case. To be precise, one needs to compute the allowed values of the mixing angle 
relevant to the solar sector the four neutrino radiative model. 

Within the standard three neutrino set up, various extensions of the Zee model have 
been presented in literature to evade the compulsion of sin^2^^Q ^ 1. Most of these include 
additional couplings and / or fields and then they can avoid bimaximal mixing and incor- 
porate the LMA solution. Here we consider the extension of the Zee model by an additional 
sterile species [|15]. This would enable us to include solutions corresponding to the LSND 
anomaly. Using the property of vanishing diagonal elements, we are able to make a pre- 
diction on the solar neutrino mixing angle involving all the three squared mass differences, 
which is compatible with the present data. We show that the model is suitable for explaining 
all the three neutrino anomalies. The precise value of 6q depends on the ratio between the 
square of the atmospheric neutrino mass squared difference and the product of the solar 
and LSND mass squared differences. Present experimental errors on these measured mass 
squared differences allow the solar neutrino oscillations to take place with MSW conversion 
either with small or large mixing. However, we find it difficult to reconcile the LOW-QVO 
solutions with this scenario, though they might not be completely ruled out. 

The rest of the paper is organized as follows. We show in section II how the near- 
maximality constraint on 6q arises from the vanishing diagonal elements of the 3x3 Zee 
mass matrix, establishing a method that can be extended to the 4x4 case. This extension 
is done in section III where we derive the result for sin^2^0 mentioned in the abstract. 
In section IV, we comment on textures of the mass matrix within this model which are 
phenomenologically viable. We then show that the active-sterile admixture in solar neutrinos 
and the LSND mass scale within these models have a common origin within the mass matrix. 
The final section, section V summarizes our conclusions. 



II. RESTRICTION ON SOLAR NEUTRINO MIXING FROM THE THREE NEU- 
TRINO ZEE MASS MATRIX 

As already mentioned, the three neutrino Zee model has a serious problem on account of 
the current neutrino oscillation data, specifically with the data from the solar sector. The 
reason for this can be seen, without going in to the details of the model, from one feature 
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of it - namely, its vanishing diagonal elements. The Zee mass matrix in three generations is 
given, in the Ue, z/^, flavor basis, by 



^0 a h\ 
a c 

c Oy 



(8) 



The elements a, b, c are three^ real parameters determining two mass squared differences 
and three mixing angles relevant to solar and atmospheric neutrino oscillation phenomenol- 
ogy. Evidently, this leads to a constrained pattern of neutrino masses and mixing which 
has to face the challenge of the emerging data. Three relations emerge from the vanishing 
diagonal terms^ of A^^-^^ if we write M^^^ = Omns diag{m,,^, m^^, m„^} Ojj^^c, : 



0, 
0, 
0. 



Because of the orthogonality of OmnS; eqs. (^), (|TD|), (|TTp lead to 



(9) 
(10) 

(11) 



m,.! + m,,, + m 



•U3 



01, - 01, 



0% 



(12) 
(13) 

Eq.(|l2l) is a trivial consequence of the tracelessness of A^^^ while eq.(|l3D implies that one of 
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the mixing angles is not independent. If Am^ = |m^ — m^. | and is positive by definition, let 



^ For three neutrinos, the Zee mass matrix has only three parameters. Thus all phases can be rotated 
away from the mixing matrix, by absorbing them in the three neutrino fields. Hence one can choose the 

orthogonal matrix Omns for mixing instead of the unitary Maki-Nakagawa-Sakata matrix Umns- 
^ Two loop corrections can induce very small nonzero diagonal entries in the Zee mass matrix of eq. 

Strictly speaking, they would change eqs. (|l^), and (|l^). However estimates of these tiny 



changes suggest that they would not significantly affect the statement of eq.(21) on the solar neutrino 
mixing angle 6*© Q. 
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us choose Am32 (= Am|^^) and Am^^ (= Am|) as the two independent mass squared dif- 
ferences. We know from condition (|T2|) that one the eigenvalues of Ai^^^ is not independent. 
The latter is chosen to be and can be eliminated. Thus one can now write 



m 



U2 



m 



(15) 
(16) 



We emphasize here that the quantities appearing in the RHS of eqs. (|T5]) and (p^, namely 
AitIq, Am^^^ are positive. These equations yield m^^^ and m^^ terms of Am^^^ and Am|. 
Let us define 



a = ^ Q . (17) 

In order to determine the solar mixing angle, we can take 623 = 6 Mm, i-e the atmospheric 
mixing angle, to be exactly maximal. Furthermore, on account of the CHOOZ p2[ exper- 
imental result, we can put Ogs = e ~ 0.1, keeping only (9(e) terms in Omns- Choosing the 
mixing angle 6^12 = 9q for solar neutrino oscillations, Omns now gets reduced to 



cos 6q sin 9q ^ 

Omns = (— e cos 6q — sin Oq)/^/! (— e sin 6q + cos 6q)/V2 1 

y (— e cos 6q + sin 6q)/\^ (— e sin 6q — cos 6q)/V2 l/V^ j 

The substitution of eq.(|T^) into eq.(|I^) leads after some algebra to 



(1 - m^Jm^^y ^ 



The physical solutions^ of m^-^, m^^ of eqs. (|^) and (|T6D yield 



Oie') . (18) 



(19) 



m 



1^1 



m 



U2 



y/1 — a + a'^ 



a 



^ We choose only those solutions which yield positive values for m^^, m^^ and sin^26'0 
are four sets with m^j ~ ±{Amj^^^/3)i {1 + 2a + 2^1" 
rriyj = ±(Am^^^/3)i (1 - 2a + 2\/l - a + a^) and m^^ 



a^) and m,. 



(20) 



There 

(-a + \/l + a + a^) vn^^, 



-{a + ^/l — a + a^) rrij^j, expressed as func- 



tions of a while m^^ is always — m^j — mi^^. The solutions are physical only for \m,. 



'Atn 



SO that y/ Am 

~ b « 0(Va;^^) and c « OiAml/./K^ 



and A/ArriQ ~ jm^j — m^J. It is not difficult to see that this implies a 



Atm/' 
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Eqs. (JT^) and (pOD lead to the following constraint on the solar mixing angle^ 







(21) 



matching with the result of Ref. |TD. Since [|, |] a = (Am|/Am^t^) ^ 6 x 10"^ and 



10 ^, we see that sin^2^0 is forced to be very close to unity {i.e | m 



i^i I 



m 



V2\ 



a situation 



disfavored by the data , though not completely ruled out. A point to note is that eq.(piD 
follows from eq.([T8|) by use of only the elements in the first row of Omns- This derivation 
does not use the fact that 62^ = ^Atm is ~ 7r/4. It is, of course, true that the extra input of 
the latter immediately forces near bimaximal mixing in the Zee model, as mentioned in the 
Introduction. Another point to keep in mind is that the texture |c| ^ |a| ~ |b|, forced by 
the neutrino mass solutions given in footnote 4, implies an approximate global Le — — L,- 
symmetry W7 . 



III. SOLAR NEUTRINO MIXING IN THE RADIATIVE MASS MODEL WITH 
FOUR NEUTRINOS 



The Majorana mass matrix with three active and one sterile species can be generally 
represented in the {z/e, z/^, u^, I's} flavor basis as 



(22) 



'^O a b 
a c e 
b c f 
d e f 

The elements a, b, c, d, e, f appearing in eq. (p2D get generated radiatively at one-loop and 
are given by 



V 



m 



(23) 



^ We should also mention that eq. (|2^) was initially derived in Rcf. ||2^ within a seesaw model which has a 
mass texture different from that of the Zee model. We thank K. R. Balaji for pointing this out. 

^ Some detailed numerical analyses of the Zee model, confronting the present data, have argued that the 
possibility sin^200 « 1 is not completely ruled out [e3, Ea, if one insists on 99% C.L. limits. 
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b = fe.K-m^)(^)F(mJ^,mJj, (24) 

c = VK-mJ)(^)FK^,mJJ, (25) 

d = (fe,f;m, + fe^f>^)/i'u F(mJ^,m2j, (26) 

e = (Vf;m, + f^efenie) yu'u F(mJ^,mJJ, (27) 

f = (M>m + fref>e) /i'u F(m^^, m^J, (28) 

where f , f (yU, /i') are dimensionless (dimensional) couplings in the model, vi 2 are two Higgs 



VEVs and m^^j, m^j are the three Higgs scalar masses, as explained in Ref.[|I^. Moreover, 
the function F(mi, 1112) is defined as 

F(mi, m2) = ^} ^ / , In ( ^ ) . (29) 
Idtt-^ — m2 \m2/ 

For simplicity, we assume the matrix elements of eq. ( p2D to be real and ignore the presence 
of any CP-violation in the neutrino sector. 

In the same manner, as shown in the previous section, one can arrive'' at the following 
relations by equating Ai^^'' with U diag{mi,j, m^j, m^^^ m^^} U""": 



m^^U^i + m^2Ue2 + m^3U^3 + m^.U^^ = 0, (30) 

m^iU^^ + m^jUjs + rn^,^l^ + m^,Uj4 = 0, (31) 

m,,U^i + m,,U^2 + ^.30^3 + m,4U^4 = 0, (32) 

m^^U^^i + m^2Us2 + m^3U^3 + m^A^l^ = 0, (33) 

where Uij are the elements of the matrix U defined in eq.(|I|). Since we ignore the phases 
present in the elements of U, we can utilize the unitarity of U to obtain the following relations 
from eq. (|0|) - (||). 

rrii^j + + m^j + vr\^^ = 0, (34) 



^ As in the three neutrino radiative Zee model, we expect that the effect of the two loop corrections to 
the four neutrino radiative model mass matrix would be small and thus of not much consequence to the 
results presented hereafter. 



10 



"el "e4 "e2 ^e4 /'Q«^ 



1)2 - 1)2 1)2 _ I |2 

'-'el '-'e4 _ '-'e2 '-'e 

U^3 - U^4 " - 

(U^2 - Ue^4)(Uf.3 - U^4) - (U^3 " Ue^4) ^2 ' ^U) 

(U^i - U^,) (U23 - U^J - (0^3 - U2,)(U2, - U^J 

(U^2 - U^4)(Ue^3 - Ue^4) " (U^2 ' Ue^4)(U^3 " ^l.) 



;U?i - U?4)(U23 - U^4) - (U^i - U2,)(U?3 - U 



(37) 



These are the four neutrino versions of eqs. (|T^), ( [I3| ) and (p!^. 

Eq. ( p^ shows us that the 2 + 2 mass pattern depicted in Fig.l can be easily made com- 
patible with this model if |m,^J ~ Irrii^J ^ |m,^3| ~ |mi,^| or if |mi^J ^ Irrii^J ^ |m,^3| ~ |miyj. 
In either case, there are three independent mass squared differences which we choose to 
be Arriji = ^it^©) ^it^43 = ^iT^Atm Am32 = ArriLSMQ. We can now perform an extended 
version of the same eigenvalue analysis as done in Section II for the three neutrino case. 



The three independent eigenvalues are chosen here to be rrii^j, mi,^, rrii^j. Eqs. (|^) and ( piq) 
now extend to 



2m^3(m^j + mj^JI = Aru/^^^, (38) 

= ^<5ND, (39) 



V2 I/i I 



Am|, (40) 



where the RHS quantities of the above equations, namely ArriLSND! Am^, Avn\^^, are positive 
as in eqs. (|T5|) and (|T6]). The physical mass eigenvalue solutions to the above three equations 
form four sets as detailed in the Appendix. From these we see that 



m,i Mf 

4 Am| Am^LSND + (Aml.J^ 



4 Am^s^o Am| ± (Aml.J^ 



+ Oia), (41) 



where are defined in the Appendix and the upper (lower) sign corresponds to sets 1, 2 
(3, 4). 
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In order to arrive at the prediction on the solar neutrino mixing angle in the four neutrino 
scenario, we can choose^ a parameterization of the 4x4 generalization of Umns matrix as 



U 



cosOq sin 6*0 e e 



V 



U^i 


U^2 


U;.3 


u 


U.i 






u 


Usl 


Us2 


Us3 


u 



s4 



(42) 



In eg. (^2]) , the solar neutrino oscillations have been chosen to be between the physical states 
|z/i > and |z/2 > where the BUGEY constraint eq.(^ has been taken care of by choosing 
Ue3 = e and Ue4 = e', both^ being - O(10~^). We keep only C(e , e') terms in the above 
mixing matrix. The rest of the parameters, while otherwise arbitrary, will be required to 
satisfy unitarity and other experimental constraints later. 

One can now determine the mixing angle in the solar sector in much the same manner 
as in the three neutrino case. Eq.® is still vahd to 0{e^,e' ^) and can be rewritten as 



Eqs. (|TI) and (g3|) together yield: 



+ 0{e\ e 2) 



(43) 



where p is defined as 



4 Am| Am^s^^j 



Eq.(^4D leads us to the following conclusions. 
• Because of eq.f^if), we have the condition 



^2 



~ J. , 



4 Am| AmfsND 



(44) 
(45) 



(46) 



* The matrix U has six independent angles and two independent phases in general. But one has the freedom 
to choose a parameterization in which the phases do not appear in Uei and Ue2- In any event, we are 
ignoring the phases. 

^ It should be noted that in the analysis presented in |7| both Ues and Ue4 are taken to be zero. 
Since mi^j/m^j ~ tan^6Q + 0{e^ , e it follows that the RHS of eq.(^l|) must be real and positive modulo 
e^, e ^ terms. 
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FIG. 2: The range of p = 2 2 ) shown on a logarithmic scale, for different solar neutrino 



4 Am^Q An.^s^Q 



solutions. 



upto a, e^, e' ^ terms. In Fig. 2 we plot the ranges of p for various proposed so- 
lutions of solar neutrino oscillations. We have chosen AiUlsmd ~ (0.2 — 2)eV^ 
1^, Am^^^ ^ (1 — 8) X 10~^eV^ [|] and, for the various solar neutrino solutions, 
Am|(SMA) ^ (2 X 10"^ - 2 x IQ-^) eV^, Am|(LMA) ^ (4 x 10"^ - 6 x lO"^) eV^ 
Am| (LOW - QVO) ^ (5 x 10"^ - 5 x 10"^^) eV^ g. From the figure, we see that 
eq. (|46D is clearly valid for large domains of SMA and LMA solutions. However, in the 



case of LOW-QVO solutions, only a small domain obeys eq.(46) making it difficult 



to be accommodated within the 4x4 radiative |]15| model of neutrino masses. A 



detailed numerical analysis of all oscillation data within the four neutrino model 
would be able to show the confidence level at which the LOW-QVO solution is valid 
given the mass matrix of eq. 



As in section III, we have expressed the solar neutrino mixing angle, 9q in terms of 
ratios of mass squared differences, though the form is different. However, unlike in 
the three generation case, the solar neutrino mixing is not forced to he very nearly 
maximal. Instead, its deviation from maximality is controUed^^ by the product of 



The special case of eq.(^ for 9q ~ 0, namely 2w ArriLSND^'^o ~ ^'^Atm discovered in Ref. ||15[ 
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the two factors in the RHS of (^41) and specifically by the LSND mass scale. If we 
take ArriLSND ~^ oo, we are back to the three neutrino result of eq.(plD. Moreover, 
the range of the ratio (Am^^^ )74Am|A 

itilsmd with present experimental errors is 
such (Fig 2) that both the small (SMA) and large mixing (LMA) MSW solutions are 
equally allowed. This is quite unlike in the three generation case. 



• The result, stated in eq. P^) , has been derived assuming the absence of any phases in 
Ail^^ and hence in U. The presence of sizeable phases {i.e of a significant amount of 
CP- violation in the neutrino sector) can change the result. 

We infer that, in so far as the prediction for the solar neutrino mixing angle is concerned, the 
radiative model with four neutrinos is on a significantly different footing as compared to 
the corresponding three generation version. Indeed, the former is better able to tackle the 
present data, especially from SNO. We discuss in the next section some phenomenologically 
allowed textures of Ai^f^ from the current data. 



IV. TEXTURES OF THE 4x4 RADIATIVE MODEL MASS MATRIX 



Our derivation of the phenomenological relation eq. (^Gf ) for the mixing angle, 6q in 
Sec. Ill made use of two inputs : (1) the 2 + 2 mass pattern implemented on the neutrino 
mass matrix of eq. ( ^21) with vanishing diagonal coefficients and (2) the form of eq. (^21) 



for the 4x4 unitary mixing matrix U. In this section, we estimate the orders of magni- 
tude of the nondiagonal entries of We do so by relating them to the mass scales 
(Aml)-'-/^, (Am^^^)-'-/^ and {^f^isNoY^^ ^he elements of U and then constraining the 
latter from phenomenological inputs. We also favor or disfavor certain mass textures of 

Let us parameterize the four real neutrino mass eigenvalues, in a way consistent with eq. 
(0), as 

nij,, = mi , = m2 , m^j = M - mi , m^^ = -M - m2. (47) 
In the 2 -\- 2 pattern, M represents the mass scale of the heavier (3, 4) pair, being 



C(y AmLSMo), while |mi^2| stand for the mass magnitudes of the lighter (1, 2) pair. Because 
of the constraint of eq. (^), we can treat (Am|)-^/^ ~ |m^ — ml\ and (AmLSNo)^ ~ |M| (see 
Fig. 1) as the two independent controlling mass scales in Ail^\ We have thus found in 
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eg. (^71) a parameterization that is convenient for the description of all neutrino oscillation 
data, obeys the tracelessness condition eq.(|3^) and follows the mass pattern depicted in Fig. 
1. 

In order to estimate the orders of magnitude of those among the nondiago- 
nal elements of Ai^^^ which are nonvanishing, we resort to the equation = 
U diag.jrrii, m2, M — mi, — M — U^. We have already put Ues ~ e, Ue4 ~ e' in eq.(E2p 
with the expectation that |e|, |e'| ~ (9(10^^) in view of eq.(@). We can also put 
U^i ^ 6, \Jfj_2 ~ S' and expect in the light of eq.(|D that \6\, \6'\ ~ 0(10"-^); indeed, the best 
fit requires them to be close to zero, cf. eq.(0). Furthermore, we can put U^3 ~ U^4 ^1/2 
on account of the observed maximal mixing in the atmospheric neutrino sector. The various 
entries of A4l^^ can now be related to the mass eigenvalues and hence to mi 2 and M as 
follows: 



a ^ mi (5 Uei - e U^s) + n^2 {6' Ue2 - e' U^4) + M (e U^s - e' U^4), (48) 

b ^ mi (Uei U^i - e U,3) + m2 (Ue2 U,2 - e' U,4) + M (e U^s - e' U,4), (49) 

c ^ mi (5 U.i - U^3 U.s) + m2 {6' U.2 - U^4 U.4) + M (U^s U.s - U^4 U.4), (50) 

d ^ mi (Uei Usi - e U,3) + m2 (Ue2 U,2 - e' U,4) + M (e U,3 - e' U34), (51) 

e ^ mi (5 Usi - U^3 Usa) + m2 {5' Us2 - U^4 Us4) + M (U^s U,3 - U^4 U34), (52) 

f = mi (Usi U.i - U,3 U.3) + m2 (Us2 U.2 - ^s4 U.4) + M (U^a U.3 - Us4 U.4). (53) 

We see from the above that in the expressions for the a, b and d, the explicit occurrence 
of e and e' in the coefficient of M pulls it down significantly; thus they are much smaller 
in magnitude compared to M. In contrast, the coefficients of M in the expressions for the 
entries c, e and f do not manifestly involve such small parameters. This is a consequence 
of the fact that the coefficients of M in all the entries depend only on the last two columns 
of the mixing matrix U eg . (^21) . We see from eqs. (|30|) - (|33|) , as also from eq. (|47|) that, in 
the limit when mi^2 0, Ug3 = Ug4, U^3 = U^4, U^3 = U^4 and = U54. In such a limit, 
when M remains as the only controlling mass scale. 



|e| ^ |e'|, (54) 
lU^sl ~ |U^4|, (55) 
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lUrsI 


~ |U.4| 


|Us3| 


~ |Us4|. 



(56) 
(57) 

We restate that we have assumed no phases to be present in our mixing matrix. In 
the hmit of neglecting terms 0{6^,5' we then have^^ from the unitarity condition 
1 = + UL + + U2 and eq.m) that 



|U;.3| = |U^4| = 1/V2. (58) 

In addition to satisfying the above conditions in the limit mi 2 0, the parameters Uss, 
Us4, U^3, U^4 also have to satisfy the off-diagonal unitarity condition: 

U^3 U^4 + U,3 U,4 + Us3 Us4 = + 0{ee'). (59) 

From eq.(p8D we see that this would mean 

Us3 U34 + U.3 U.4 ~ (60) 

Eqs. ([55|), (0) and { ^7\ ) admit four possibilities corresponding to the signs each of the mixing 
matrix elements can take. First, we rewrite the condition (^) in terms of two possibilities as 
U/^3 = U^4 or U^3 = — U^4. The remaining possibilities form sub-cases of these two cases. 
Each of these (sub) cases is further constrained by the off-diagonal unitarity condition of 
eg. (|5PD . After some algebra, we realize that three possible cases are consistent with the 
unitarity conditions. Each of them leads to a different texture which we detail below. In all 
the following cases, we have a, b, c ~ 0(mi 2). 

• Case A Here, we estimate the following orders of magnitude for the other entries of 

c fti 2 M U^3 U,3 = M'; e ^ 0(mi,2); f ^ 2 M (J^s U.3 = M". (61) 

This case comes either of the choices: U^3 = U^4, Ut-3 = — Ur4, Us3 = Us4 or U^3 
= — U^4, U^3 = U^4, Us3 = — Us4. The mass texture would be of the form : 



12 



It follows from eq.(p^) that, when vni 2 ^ 0, i^s and 1^4 form a Dirac neutrino. So it is not surprising 
that we get maximal mixing in the 3, 4 sector in this limit. 
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^ a b d ^ 

a M' e 

b M' M" 

d e M" 

The active sterile admixture A is now given by 

f2 



V 



(62) 



A ^ 1 



-2 ■ 



(63) 



We see that^^ 2 U53 < 1. Thus we can parameterize f/c as cos (3 in which case 
A = sin^ /3. 

Case B Here the various magnitudes that c, e, f can take are given by 

c ^ 0(mi,2); e ^ 2 M U^3 Uss = M'"; f ^ 2 M Us3 U.s = M". (64) 

This situation can arise with either of the choices : U^3 = U^4, U^s = 11^4, Uss = — Us4 
or U^3 = — U^4, Ur3 = — Ut-4, Us3 = Us4. The mass matrix would now take the form: 






a 


b 


d 


a 





c 


M'" 


b 


c 





M" 


d 


M'" 


M" 






(65) 



The active sterile admixture in this case is given by 



(66) 



Once again, f^/e^ = 2 and can be parameterized as cos^(3 so that A = sin^ (3. 
Case C Here we have 

c ^ 2 M U^3 U.3 = M'; e ^ 2 M U^3 ^ss = M'"; f ^ 0(mi,2), (67) 
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The unitarity condition 1 = Uji + U^j + U^j + U^4 and eq. (|57|) imply that f^c^ ~ |Us3| ~ 1/2. 
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arising from either of the choices : U^3 = U^4, 11^3 = 11^4, Uss = Us4 or U^js = — U^4, 
Utau3 = — Ur4, Us3 = — Us4. In this casc the mass matrix has the following texture 



a b d 

a M' M'" 

b M' f 

d M'" f 







(68) 



The active sterile admixture is now given by 



, (69) 
c2 + e2 ^ ^ 

The form A = sin^ (3 can now be achieved by parameterizing e/c = Urs/Uss as tan f]. 



We note that this texture has been recently been studied in Ref. ||2^ on the basis of two 
pseudo-Dirac neutrinos and an approximate global Le + — — Lg symmetry. 
A slight variation of eq. (|6^) with f = and a nonzero fourth diagonal element has 



been proposed in Ref. on the basis of a paired pseudo-Dirac structure and the same 
global symmetry, but not within a Zee-type radiative model. The textures of our cases 
A and B strongly violate the global Le + — L,- — Ls symmetry. 

All the above textures give rise to the 2 -\- 2 pattern along with an admixture which 
in general is between zero or one. The choice Uss = or U53 = 1/2 will yield the extreme 
results A = or 1 respectively, but such a choice is not favored by the experimental data [0. 
In all the above three cases the LSND scale appears in two pairs of elements in Ail^'^ if (3 
is nonzero. This is not surprising since, for a nonvanishing admixture A, the sterile species 
Us and the tau neutrino z/^ appear in both the lighter (1,2) and the heavier (3.4) pair of 
physical states. 

V. CONCLUSIONS 



In this paper we have studied the constraints on the neutrino mass matrix Ai^, for a 



radiative model |15] of three active and one sterile neutrino species from all the current 
data on neutrino oscillations. We have established a method to show how the solar neutrino 
mixing angle is related to the observed squared mass differences. This leads, in the three 
neutrino case to a near-maximal 6q, not in favor with the latest data. In the four neutrino 
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Set 1 


Set 2 


Set 3 


Set 4 


m,, = N+/D+ 


m,, = -N+/D+ 


= N^/D- 


m^i = -N^/D- 


m,, = N+/D+ 


m,, = -N+/D+ 


= N^/D- 


= - N2 /D- 


m,3 = N+/D+ 


m.3 = -N^/D+ 


m^3 = N3"/D- 


m^3 = -N3/D- 



TABLE I: Neutrino mass solutions 



case, however, we derive the result sin^2^^0 ^ 1 — [(Am^^^)^/ (4 ArriLSND^'^o)]^ which allows 
a phenomenologically acceptable value for Oq within the allowed ranges of the three squared 
mass differences. The above result is compatible with the present solar neutrino data when 
MSW solutions with either large or small mixing are considered. However, it is difficult to 
reconcile LOW-QVO solutions within these models. 

We have shown then that the radiative model with four neutrinos offers several phe- 
nomenologically viable textures of M-^f^ ■ This is quite unlike in the three neutrino case 
where maximal mixing in the atmospheric sector forces the experimentally disfavored bi- 
maximal texture. 

Acknowledgments: We wish to thank Anjan Joshipura, Uma Sankar and Mohan Narayan 
for fruitful discussions. 

APPENDIX 

Here we present the neutrino mass solutions of eqs.(|38D - (0). There are four sets of 
solutions. These can be conveniently expressed in terms of numerators N^2,3 and common 
denominators defined as follows. 

N+ = - 4 Am^Ls^D Am^ + (Am^^,)^ - 2 Am^ Am^^^ - 3 (Am|)^ 
N+ = 4 Am^Ls^D Am| + (Am^^,)^ + 2 Am| Am^^, + (Am^)^, 
N3+ = 4 Am^Ls^D Ami,^ + (Ami,^)^ + 2 Am| Am^,^ + (Am|)2, 
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and 



D^ 



l\fl [{(Am 



Atm J 



[Am|)2} (Am^^, + 2 Am^s^^ + Am| 



= 4 Am^s^o Am^Q + (Am^^,)^ + 2 Am^ Am^^, - 3 (Am|)^ 



2 \2 



4 AmfsND Am^ + (Am^,,)^ - 2 Am^ Am^,, + (Am^) 



= 



D 



4 Am^s^D Am^^^ 



[Am2^,)2 - 2 Am| Am^^, 



(Ami) 



2n2 

0J ' 



2V2 [{(Am 



Atm J 



(Am|)2} (Am2^, + 2 Amfg^p - Am| 



Note that N|~ (i = 1,2,3) and D~ can always be obtained from the corresponding and 
by putting Am|, — Am|. The four sets of neutrino mass solutions with xWy^ always 
being —my^ — xx\y^ — xx\y^, are shown in table I. 
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